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1 Introduction 

The neutrino deep-inelastic scattering (DIS) is continuing to serve as the classical 
tool for probing the nucleon structure (for reviews, see e.g. Refs. lil* nJ). In the last 
few years some progress was made in more detailed experimental and theoretical 
studies of the behaviour of the cross-sections of vN DIS, and in the extraction of 
the structure functions (SFs) XF3 and F2. Among the data that are currently under 
active analysis are the ones provided by the CCFR/NuTeV collaboration at the 

Fermilab Tevatron (see e.g. Refs. 0/ a/), the experimental results of theJINR-IHEP 

Neutrino Detector collaboration EP collected some time ago at the IHEP (Protvino) 
U-70 proton synchrotron, and the preliminary data of the CHORUS collaboration, 

obtained recently at the CERN SPS §) 0). 

The kinematical regions currently available for cross-section measurements 
are shown on the plot of Fig. |l|, taken from Ref. El'' . This figure clearly shows that the 
above-mentioned three experiments were performed in different kinematical regions, 
which overlap in part only. Thus they provide complementary information about 
the behaviour of SFs in different regimes. Moreover, additional more precise data 
for uN DIS cross-sections can be obtained in the future at neutrino factories. If the 
energy of the neutrino beam is fixed at E^, = 50 GeV, experiments should penetrate 
into the physical region that was added to Fig. 1 in Ref. B' . This region overlaps in 
part with those where the above-mentioned three experimental collaborations were 
working. Therefore, the studies of the data on uN DIS characteristics available at 
present can be important milestones in the planning of more precise DIS experiments 
at neutrino factories i-^K 

2 Discussions of some new experimental results 

Recent interesting experimental news came from the model-independent re-extraction 
of the behaviour of the F2 uN SF □/ from the CCFR'97 data ™ . The re-analysis 
of Ref. a/, which does not affect previous CCFR'97 xF^ results, removed the widely 
discussed discrepancy that existed at x < 0.1 between the behaviour of CCFR'97 
F2 B* and that obtained by the NMC collaboration liil'' from the process oi fxN DIS. 
In addition to the new extraction of F2 from the differential cross-sections of CCFR, 
the first measurement of AxF^ = xF^ — xF^ was also performed a) . However, it is 
important to stress that none of the considered non-perturbative, perturbative, and 
theoretical effects, including charm production, are still unable to describe properly 
the experimental results obtained for AxF^ \=^ . 

The information about the corrected behaviour of CCFR F2 may be really 
important for the continuation of the work of the CHORUS collaboration. Indeed, 
in the kinematical region where both sets of uN DIS data overlap, the preliminary 
CHORUS results for xF^ and F2 \-^ agree with the ones provided by the CCFR 
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Figure 1: Kinematic regions in x-Q^ for cross-section measurements in deep inelastic 
ep scattering, v scattering and for triple differential jet cross-section measurements 
in pp collisions (from Ref \^^ ). 



collaboration in 1997 H/ . Therefore, the preliminary CHORUS F2 data should show 
a pattern identical to that found in the CCFR'97 analysis ^^ , i.e. exceeding by 10- 
15% F2 NMC measurements at x < 0.1. This excess is beyond the existing statistical 
and systematic errors of discussed DIS experiments. Moreover, the inclusion of 
these wrong CCFR'97 F2 points into the next-to-leading order (NLO) QCD fits, 

performed with the help of the DGLAP method \±^ , leads to the erroneous low-x 
behaviour of the gluon distribution xG{x, 9 GeV^) ~ x*'^ with be = 0.0092 ± 0.0073 
z3l . It is in evident contradiction with the number, obtained previously from the 
NLO combined analysis of the data from HERA and the CERN SPS LL^, namely 
be = -0.267 ± 0.043 at Qg = 9 GeV^ Taking into account new CCFR model- 
independent extractions of -F2 eI/, it seems worth while to perform more careful 
studies of the preliminary CHORUS results. Moreover, it is rather interesting to 
try to verify from the CHORUS data the experimental behaviour of AxF^, found in 
Ref. I). 

It should be stressed that the CHORUS experiment has an attractive fea- 
ture. Indeed, as can be clearly seen from Fig.l, it provides information about uN 



Nuclear shadowing ejfect for F2 and F3 




0.0001 



0.001 



X 



Figure 2: The ratios of a heavy target to the free nucleon SFs R2 
R3 = F^/F^ calculated for ^^Fe nucleus in the region of small x and Q 

(the figure of Kulagin from Ref. \11jJ ) 
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DIS SFs in the region of rather low Q^ and low x, which complement in part the one 
where the CCFR'97 data were extracted. In this kinematical domain, theoretical 
contributions of l/Q"^ and nuclear corrections can play an important role. Leaving 
for a while the discussions of power-suppressed terms, we stress that the CHORUS 
collaboration was using a lead target, while the CCFR target is made of iron. Pos- 
sibilities are therefore really open to study nuclear effects in neutrino DIS; as shown 
in calculations reported in Ref. li^J, these effects can be of great importance. A 
comparison of these effects for the cases of F2 and xF^ neutrino DIS SFs is depicted 
in Fig. ^ constructed for the detailed work of Ref. c}^ . 

Another interesting possibility of DIS experiments is the extraction, from 
their chraracteristics, of non-perturbative power-suppressed terms and the values 
of as{Mz)- This question was considered in Refs. lill/ c^) in the process of the 



next-to-next-to-leading order (NNLO) QCD fits to different data, and in Ref. cS*, 
while performing NLO fits to the experimental results for charged-leptons DIS SFs. 
The most recent outcomes of the NNLO analysis of the CCFR'97 xF^ data c^ will 
be discussed in the next section. Here it is worth while emphasizing that the works 
of Refs. LLlI/ c}3 \£^) E2I/ agree in their conclusion that the inclusion of the NNLO 
QCD corrections into the fits has a tendency to decrease the extracted values of 
non-perturbative 1/Q^-terms. Whether this is a general theoretical feature (see e.g. 

Refs. ^21/ EH/ ) or it is related to the lack of precision of the analysed data might be 
clarified in the future, if the ideas of more detailed experiments on neutrino DIS at 
neutrino factories \^ are realized. 

It should be noted that the situation at NLO is more transparent. Indeed, 
the Jacobi polynomial fits of Refs. \=-^ \=^ \=^ demonstrated that it is then pos- 



Table 1: The results of the fit, in Ref. u) , of the IRR model to the data from 
different neutrino experiments. The value of y^ over the number of points (np) is 
given. 



Experiment 


Ai GeV2 


X^/np 


IHEP-JINR 


0.69 ± 0.37 


3/12 


CCFR'97 


0.36 ±0.22 


253/222 



sible to determine both 1/Q^ twist-4 effects and as{Mz) from tlie CCFR'97 xFs 
experimental results, cut at Q^ > 5 GeV . Moreover, the NLO DGLAP fits of 
Ref. li^ confirmed this feature, using both xF^^ and a cut at x > 0.1 on F2 vN 
CCFR'97 data. A similar conclusion was made in the process of NLO DGLAP fits 
of intermediate-energy charged leptons data c3i) . However, it is obvious that it is 
simpler to detect non-perturbative effects in a smaller Q^-region. 

The most recent example of such a determination was given by the NLO 
DGLAP analysis 13/ of DIS neutrino data from JINR-IHEP Neutrino Detector col- 
laboration, which collected rather large statistics (5987 neutrino and 741 antineu- 
trino charged-current events) in a rather low-energy region 0.55 GeV^ < Q^ < 
20 GeV^. The non-perturbative l/Q^ correction to the perturbation theory (PT) 
behaviour of xF-^ was parametrized using the infrared renormalon (IRR) model of 

Ref. E3), namely 



xF^{x,Q')=xFp{x,Q') + 



h{x) 



where 



h{x) = A 



^,'^C^^^^(z)p^six/z,Q') 

X Z 



(2) 



(jHiR ^g |-]^g IRR model coefficient function, xpns{x,Q'q) = Ax^'^^{l — xY'^^ is the 
boundary parton distribution, defined at Ql = 0.5 GeV^ at NLO, and A2 is the 
normalization parameter, which should be determined from the fits. It is sometimes 
expressed through the parameter A3 as 



A' 



2CV . 2 



(3) 



where Cp = 4/3 and /?o = (H ^ 2/3/) is the first coefficient of the QCD /3-function. 
Fixing as{Mz) = 0.118, which corresponds to its world-average value, A3 was ex- 
tracted from the fits to the IHEP-JINR Neutrino Dectector and CCFR'97 data D. 
The results are presented in Table 1, taken from Ref. W . Note that the two expres- 
sions for A3 from Table 1 are comparable within the errors. Averaging the numbers 
for A3 and transforming them to A'2, the authors of Ref. ^ obtained the following 
value: 



A'2 = -0.130 ± 0.056 (exp) GeV^ 



(4) 



where the error includes both statistical and systematic experimental uncertainties. 
It is in agreement with the value extracted from the NLO Jacobi polynomial analysis 
of the CCFR'97 xF^ behavioun^cut at Q^ > 5 GeV^ 113) HI HI). Indeed, at NLO, 
the most detailed fits of Ref. c^) give : 

^2 = -0.125 ±0.053 (stat) GeV^ . (5) 

Note, however, that the error in Eq.(5) does not include the systematic uncertainties. 
Therefore, Eq. (4) is the most precise up-to-date value of the IRR model parameter 

3 New QCD fits to CCFR xFg data: NNLO and beyond 

We can start the discussion on the phenomenological application of some new N^LO 
perturbative QCD results on the coefficient functions of odd moments of xF^ and 
on the NNLO approximations for the related anomalous dimensions ^21/ (which are 
complementary to those obtained in Ref. ^MJ in the case of even moments for the 
F2 SF of charged- leptons DIS), in combination with the NNLO expressions for the 
coefficient functions Eiil' that were recently confirmed in Ref. ^-^> . 

3.1 The application of the Jacobi polynomial method 

It is appropriate, at this point to recall the basic ideas of the Jacobi polynomial 
method cl^ which was developed in Refs. Etf EJ/ and was previously used in the 



analysis of the BCDMS charged-leptons DIS data at NLO ^'', and in the non-singlet 
approximation at NNLO ^21/ . In the case of the analysis of the CCFR xF^ data, the 
Jacobi polynomial method was applied at NLO in Refs. cJV l22l/ and proved useful 
for performing fits at the NNLO 113) I19D E^D and approximate N^LO levels, H^) 1^) 
with and without twist-4 corrections (see discussion below). 

This method allows the reconstruction of the SF (say xF^) from the finite 
number of Mellin moments, namely 



xF3^-(x,g^) = x"(i-xf Ee^^(x)Ecf(«,/?)M™53(g^) + ^ , (6) 



E0n^(^)E 
n=0 jr=0 



where 6^'^ are the Jacobi orthogonal polynomials with parameters a, /3; c^ (a, f3) is 
the combination of Euler F-functions. It increases factorially with increasing n. The 
Mellin moments M™^ include information on the 1/Q^ target mass corrections and 
are defined as 

<r (Q') = Jl ^-'Fr{x, Q-)dx + "^^^^M^^^^fM) ■ (7) 



Table 2: The NNLO results 
mined, in Ref. ^cn> and their 
new ones are marked by hold 



of the parameters A, b, c of the model 
comparison with the values obtained in 
type. 



for xF^ deter- 
Ref E3) . The 



Oidei/Nmax 


Ql 


A 


b 


c 


xVnp 


NNLO/6 


5GeV2 


4.25±0.38 


0.66±0.03 


3.56±0.07 


78.4/86 


NNLO/9 




3.73±0.68 


0.63±0.05 


3.52±0.08 


72.4/86 


NNLO /6 


10 GeV^ 


4.50±0.36 


0.65±0.03 


3.73±0.07 


76.3/86 


NNLO/9 




4.21±0.35 


0.63±0.03 


3.73±0.07 


74.2/86 


NNLO/6 


20 GeV^ 


4.70±0.34 


0.65±0.03 


3.88±0.08 


77.0/86 


NNLO/9 




4.49±0.25 


0.63±0.02 


3.89±0.06 


75.8/86 


NNLO/6 


100 GeV^ 


4.91±0.28 


0.63±0.02 


4.11±0.10 


80.0/86 


NNLO/9 




4.74±0.32 


0.61±0.02 


4.14±0.09 


77.8/86 



The contribution of the twist-4 terms to Eq.(6) is parametrized with the help of the 
function h{x). It will be neglected for our first stage of discussions. 
Fixing now the behaviour xF^ at the initial scale Qq as 

xFrix,Ql)=AiQl)x'^'^o)(^l-:,y(Ql)(^l + ^(^Ql)) , (8) 

calculating the related Mellin moments and transforming them to experimentally 
accessible regions with the help of the renormalization group technique at LO, 
NLO, NNLO and approximate N^LO (the explicit formulae for the renormalization 



group evolution can be found in Ref. I£2I/), substituting the renormalization-group- 
improved expression for M™^{Q'^) into Eq. (6), and performing the fits to the 
experimental data, it is possible to determine 5 parameters, namely A, b, c, 7 and 
AijL; this enters the QCD coupling constant as, defined up to N'^LO, by means of 
the solution of the renormalization group equation for the explicitly known 4-loop 



approximation of the QCD /^-function cM) . In Table 2 the new NNLO results of 
Ref. \=^ for the parameters A, b, c of the model of Eq. (8) are presented. One can 

notice that, although these results agree with those of Ref. cl3 within the statistical 
errors, the central values of the new numbers are over 0.03 lower, and the new fits 
have smaller values of x^- 

It should be stressed that the new multiloop calculations of Ref. 1^'' allow 
us to use more moments in Eq. (6), namely n = 13, which corresponds to fixing 



iV„ 



9. All this information was effectively used in Ref. 



Note that the 



previous, similar xF^ fits of Refs. li-i' \^) were made in the case of ra = 10 and Nmax = 
6; they were using the approximate NNLO expressions for non-singlet anomalous 
dimensions, obtained from exact NNLO expressions for the anomalous dimensions 
of even non-singlet moments of F2, calculated in Ref. 1^''. Thus the considerations 
of Ref. c£P contain less uncertainties than the previous analysis of Refs. cJl) \±3) . 
The NNLO and N^LO results for A^^, obtained in Ref. i=^l in the process of twist-4 



MS^ 



Table 3: 



The Qo o^^'^d N„ 



dependence of AjJ^ ( in MeV) from Ref. L~d/ . The values 



of X o,f^ presented in parenthesis. 



MS 



N 


Q'o (GeV^) 


5 


8 


10 


20 


50 


100 


6 


NNLO 


297±30 


314±34 


320±34 


327±36 


327±35 


326±35 






(77.9) 


(76.3) 


(76.2) 


(76.9) 


(78.5) 


(79.5) 


7 


NNLO 


326±34 


327±35 


327±35 


326±36 


327±36 


328±35 






(75.9) 


(76.7) 


(77.1) 


(78.1) 


(78.8) 


(78.7) 


8 


NNLO 


334±35 


334±35 


333±35 


331±35 


328±35 


328±35 






(74.3) 


(75.7) 


(76.2) 


(77.4) 


(78.3) 


(78.5) 


9 


NNLO 


330±33 


332±35 


333±34 


331±37 


330±35 


329±35 






(72.4) 


(73.6) 


(74.7) 


(75.8) 


(76.7) 


(77.8) 


6 


N^LO 


303±29 


317±31 


321±32 


325±33 


325± 33 


324±33 






(76.4) 


(75.6) 


(75.7) 


(76.6) 


(78.0) 


(78.7) 


7 


N^LO 


328±32 


326±33 


325±33 


322±33 


324± 33 


324±33 






(76.2) 


(77.0) 


(77.3) 


(78.2) 


(78.5) 


(78.2) 


8 


N^LO 


334±33 


329±33 


327±34 


324±34 


323± 34 


324±34 






(74.8) 


(76.2) 


(76.6) 


(77.4) 


(77.3) 


(77.2) 


9 


N^LO 


330±31 


329±34 


329±32 


325±33 


325± 32 


325±33 






(73.3) 


(74.6) 


(75.7) 


(76.4) 


(76.7) 


(76.8) 



independent fits, are presented in Table 3. One can see that the application of new 
information from Ref. c£P, which allowed N^ax in Eq- (6) to go from 6 to 9, leads to 
better stability of Ai^ with respect to changes of the initial scale Qq, and decreases 
the values of x^- 

It is interesting to compare the NNLO results for ^(Qo) from Table 2 , 
which are almost Qo-independent, with the calculations of the small-x asymptotic 
behaviour of non-singlet contributions to Fi and the spin-dependent SF gi performed 
in Ref. ^^ in all orders of 1-loop expression for as, using in part the approach 
developed in Ref. c^^ . In the process of calculations of Ref. ciJ'' the following 1-loop 
formula for a, was used 



a< 



is) 



Att 



An 



(3o ln(-s/A2) f3o[ln{s/A^) - tn] 



(9) 



Att 
To 



ln(s/A2 



+ 



i-n 



ln'(s/A2) + 7r2 ln'(s/A2) + tt 



where A = A^-^ ^^ = 0.1. Note that the idea of taking into account the effects of 
TT^-terms in the perturbative expansion parameter is not new. It was previously 
used in a number of works on the subject c^f" crP (for recent applications see in 

particular Ref. \^) and Ref. e3/ where other related works were discussed as well) . 
The solution of the corresponding equations, which will not be presented 
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Figure 3: Dependence of ujq on rj = ln(/i^/A^). 1: for F-l^'^ , 2: for g^^ , 3 and 4- 
for F^^ and g^^ , respectively, without accounting for the n'^-terms (from Ref. aJjJJ, 
The results of the fits of Ref. 1^ are added. 



here, and the apphcation of Eq. (9) allowed the authors of Ref. cil'' to obtain the 
following small-x asymptotic behaviour of F^^ and g^^ : 



^1^ 



NS 



9i 



NS 



"^0 



"^0 






4+)/2 



4-'/^ 



(10) 



where the powers cOq and cuq are drawn in Fig. 3, taken from Ref. ciP. This plot is 
supplemented with the bounds on the values of cjq = 1 —b{Ql), which come from the 
results presented in Table 2 and obtained in the process of the fits to the CCFR'97 
data for the non-singlet SF xF^, performed in Ref. caJ for Qq>5 GeV^. These 
values are lying in the area limited by straight lines, with the centre at ujq = 0.37 
shown in Fig. 3. It is worth trying to understand in more detail the possible 
relations between the results of Ref. c^/ and those of Ref. c2l/. Note, in particular, 
that the comparison seems to be more legitimate in the region rj > 8, namely for 
Q^ > fi^ = 30 GeV^, than for low values of rj, where the 7r^-effects are playing the 
dominant role. This region is not considerably affected even by the transformation 

from A^^) Ri 0.1 GeV, used in Ref. cllJ, to A^^^ ^ 0.4 GeV, which follows from the 
Qg-independent numbers of Table 3, obtained when Nmax is fixed to 9 in Eq. (6). 



Consider now some other results of the work of Ref. 1^21/ , and in particular 
the extraction of twist-4 contributions and the value of as{Mz) at various orders 
of perturbation theory. To model the 1/Q^-term h{x) in Eq. (6), three approaches 

were used in Ref. c^ . The first one is the IRR model of Ref. cJj) (see Eq. (2)). 

After taking the Mellin moments from Eq. (2) and applying the NNLO 
and N^LO fits to the CCFR'97 xF^ data, the reductions of the NLO value of A'^, 
presented in Eq. (5), were observed. At the NNLO and N^LO the expressions for 



(11) 



A2 become comparable with zero, within the statistical errors c^ , namely 



NNLO 
N^LO 



A, 

a: 



-0.013 ±0.051 GeV^ 
0.038 ±0.051 GeV^ 



However, the related as{Mz) results were determined in Ref. c^ with reasonable 



errors: 



(12) 



NLO : a,{Mz) = 0.120 ± 0.002 (stat) ± 0.005 (syst) 

± 0.002 (thresh) +°:°J° (scale) 
NNLO : a,{Mz) = 0.119 ± 0.002 (stat) ± 0.005 (syst) 

± 0.002 (thresh) +°:°°^ (scale) 
N^LO : as{Mz) = 0.119 ± 0.002 (stat) ± 0.005 (syst) 

± 0.002 (thresh) +°:°°2 (scale) 

where the first theoretical uncertainty is due to the ambiguities of taking into account 
threshold effects while transforming the results for A^ to a world with f = 5 

o MS •' 

numbers of active flavours (for a detailed explanation of how this uncertainty was 

fixed using the matching conditions of Ref. e^ , see Ref. f^) and references therein) 
and the scale-dependence uncertainty was determined by choosing the factorization 
and renormalization scales n'^p = ix\ = /^Ijo k and varying k in the conventional 
interval 1/4 < k < A. One can notice the drastic reduction of the scale-dependence 
uncertainties as a result of adding NNLO and N^LO perturbative QCD corrections 
into the fits, tabulated in the case of / = 4 in Ref. 1^'' (note that at the N^LO the 
contributions to expanded anomalous- dimension terms were modelled using [1/1] 
Fade approximants). 

The results for as{Mz), presented in Eq. (12), should be compared with the 
ones obtained from the twist-4 independent Jacobi polynomial fits to the CCFR'97 
data at Nmax = 9 ^^U , which give 



NLO : a,{Mz) 

NNLO : as{Mz) 

N^LO : as{Mz) 



0.118 ± 0.002(stat) ± 0.005(syst) 
± 0.002(thresh)l°:°°^(scale) 

0.119 ± 0.002(stat) ± 0.005(syst) 
± 0.002(thresh)l°:°°^(scale) 

0.119 ± 0.002(stat) ± 0.005(syst) 
± 0.002(thresh)l°:°°^(scale) 



(13) 



Notice that the effective minimization of the twist-4 contributions at the NNLO and 
N^LO (see Eq. (11)) is leading to rather closed NNLO and N^LO values of as{Mz), 
which were obtained from the fits with and without l/Q^ corrections. 

It is worth stressing that errors on the scale dependence of the NLO and 
NNLO results from Eq. (13) have definite support. Indeed, they are in agreement 
with the independent estimates 

Aas{Mz)NLO =t'o'ool , Aa,{Mz)NNLO=t'oZll , (14) 

obtained in Ref. BiJ/, which use the model constructed in this work for the NNLO 
NS DGLAP kernel. 

In order to study the second possibility of modelling 1/Q^-effects using the 
parametrization of h{x) by free constants hi = h{xi), where Xi are the points in 

the experimental data binning, 9 parameters hi were used in Ref. c:3) . This choice 
distinguishes new fits from the ones performed in Refs. LLd'' c3) , where 16 variables hi 
were used. The minimization of the number of free parameters was motivated by the 
works of Refs. cS li^J , where it was demonstrated that a decrease in the number of 



fitted high-twist parameters decreases the correlation between their errors and make 
their extraction more reliable (the problems of estimating theoretical uncertainties 

in the case of the choice of 16 free parameters hi were also discussed in Ref. ^''). 
The choice of a smaller number of hi results in a more reliable description of the x- 
shape of h{x) for the fits to the CCFR xF^ data. As in the process of the analogous 

fits of Refs. lill' c3) , the LO and NLO x-shapes of h{x) obtained in Ref. c£P are 

in agreement with the prediction of the IRR model of Ref. clU) . The new NNLO 



and N^LO results of Ref. cB) , in agreement with the above-discussed tendency to an 
overall minimization of the extracted contribution of h{x), reveal some new feature, 
namely an indication of an oscillating-type behaviour of h{x) around x = 0, albeit 
with rather small amplitude. 



The third model of 1/Q^-corrections, considered in Ref. ^^, is directly 
expressed in terms of Mellin moments, namely 

Mr{Q')=n^M^HQ') (15) 



with the free parameter B2. It is identical to the model used in Ref. ^-J' for fixing 
theoretical uncertainties of the extraction of as{Mz) at NNLO with the help of the 
Bernstein polynomial technique, which will be discussed below. 

In Ref. c^ it was shown that the precision of the CCFR'97 xF^ data 
allows a determination of the value of B2 together with as{Mz) both at LO and 
NLO. However, as in the cases of the previous two models for l/Q^-corrections, 
considered in Ref. c2l/ , NNLO perturbative QCD effects screen the contribution of 
non-perturbative l/Q^-corrections, defined through Eq. (15). It should be recalled. 



in these circumstances, that the QCD fits of Ref. K'' to BEBC-Gargamelle 1^ and 

CDHS ES neutrino DIS data, performed over 20 years ago, did not allow a dis- 
crimination between l/Q^ and the logarithmic description of scaling violation to be 
made. Therefore, it is possible to conclude that present neutrino DIS data now have 
become more precise. Indeed, their analysis shifted the effect of perturbative screen- 
ing of l/Q^-corrections from LO to NNLO. The next generation of more detailed 

tests of QCD in neutrino DIS is now on the agenda c-^ . 

3.2 The application of the Bernstein polynomial method 

In this part of our mini-review the basic steps of the Bernstein polynomial approach, 
proposed in Ref. 1^ and recently used in the process of NNLO fits to the CCFR'97 
xF^ data in Ref. i=^ , will be recalled. The basic constructions of this approach are 
the Bernstein averages for the xF^ SF : 

F^liQ'') = C dxpnk{x)xF;{x,Q'') , (16) 

where Pnk{x) are the Bernstein polynomials, which can be presented, when k < n, 
in the following form: 

n-k ( — lY 

Pnk{x)=p{n^k)j:^ ^^_^^\^_^^/ ^'^'^^' , (17) 



where p(n, k) is defined as (see e.g. til/): 

2(n-^!r(n + |) 

Using Eqs. (16)-(18), it is possible to express the Bernstein averages for xF^ through 
xF^ odd Mellin moments as: 

n-k f — iV 

f:^{Q') = Pin. k) Y: ^,^^^_^^_^^, M(,,+2m),F3(g^) • (19) 

The next step is similar to the one used in the process of applications of 
the Jacobi polynomial technique. At some initial scale Qq, xF^ can be parametrized 
through Eq.(8), several odd Mellin moments of xF^ at Q^ defined and then trans- 
formed using NNLO renormalization group equations at the appropriate values of 
Q^, which enter into the kinematical region of the analysed experimental data. Form- 
ing now Bernstein averages of Eq. (18) it is possible to fit them to their experimental 
values. In the case of the CCFR'97 xF^ data, fits were made in the kinematical re- 
gion 7.9 GeV^ < Q'^ < 125.9 GeV^ 0). The following numbers for a^(Mz) were 
obtained IriP from these fit 



s 



NLO : as{Mz) = 0.116 ± 0.004 (exp) (20) 

NNLO : as{Mz) = 0.1153 ± 0.004 (exp) , 



The final NNLO expression, which includes the estimates of some theoretical uncer- 



tamties, is 



21) 



NNLO : a,(Mz) = 0.1153 ±0.0041 (exp)± 0.0061 (theor) , (21) 

It is worth while to mention that, despite the qualitative agreement, the central 
NLO values of Eq.(20), obtained with the help of the Bernstein polynomial tech- 
nique, are lower than the existing determinations of as{Mz) from the CCFR'97 
xF^ data, which result from the NLO DGLAP analysis B* liS and the application 
of the Jacobi polynomial technique \Lll) ]±A) c^ . Moreover, at NNLO, the result 
of Eq. (21) intersects with the NNLO determination of as{Mz) of Ref. ^^ (see 
Eqs.(12) and (13)) within existing errors only. The comparison between the results 
of the Jacobi and Bernstein polynomial determinations of as{Mz) and of the related 

theoretical uncertainties was presented in Ref. \=^ . In the process of these studies, 
definite disagreements were revealed between some results of the works of Ref. 1^'' 
and Ref. c^) . The origin of these disagreements is unclear at present and stimu- 
lates a more detailed analysis of the NNLO realizations of the Jacobi and Bernstein 
polynomial approaches. Note, however, that the definite choice of the scale param- 
eter in the Jacobi polynomial fits leads to improving the agreement of the results of 



applications of the two methods c3' . In view of this observation, it is possible that 
the results of Ref. ^il' contain larger theoretical uncertainties due to the neglect of 
scale-dependence ambiguities. On the other hand, contrary to the Bernstein poly- 
nomial analysis, the NNLO Jacobi polynomial fits of Ref. ^) also used approximate 
information about the values of the NNLO corrections to anomalous dimensions of 
even moments of xF^. It should be stressed that this approximation can be elimi- 
nated after completing the program of explicit calculations of NNLO contributions 
to non-singlet DGLAP kernels, which is now in progress ^K As to the current 
applications of the DGLAP method in the concrete NNLO fits to DIS data, they 
can in principle be based on the machinery of the Bayesian treatment of systematic 



errors of DIS data (see e.g. Ref. Ii3/) and the approximate NNLO models of DGLAP 
kernels, constructed in Refs. 1^'' 
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